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Abstract 

This paper gives a survey of a diagrammatic approach for fermionic 
pairing fluctuations, which are relevant to the BCS-BEC crossover real- 
ized with ultracold Fermi gases. Emphasis will be given to the physical 
intuition about the relevant physical processes that can be associated 
with this approach. Specific results will be presented for thermody- 
namic and dynamical quantities, where a critical comparison with al- 
ternative diagrammatic approaches will also be attempted. 



[To be published as a chapter in BCS-BEC Crossover and the Uni- 
tary Fermi Gas (Lecture Notes in Physics), edited by Wilhelm Zwerger 
(Springer, 2011).] 



1 Introduction 



The BCS-BEC crossover has been of considerably interest over the last sev- 
eral years, especially after its experimental realization with ultracold Fermi 
( 6 Li and 40 K) gases pQ. By this approach, a continuous evolution is sought 
from a BCS-like situation whereby Cooper pairs are highly overlapping, 
to a BEC-like situation where composite bosons form out of fermion pairs 
and condense at sufficiently low temperature. Here, reference to composite 
bosons stems from the fact that the temperatures of formation and conden- 
sation are in this case comparable with each other, in contrast with more 
conventional point-like bosons for which the two temperatures are quite dif- 
ferent (reflecting the fact that their internal structure has no relevance to 
problems related to condensation). Accordingly, a theoretical description 
of composite bosons should take into account not only their overall bosonic 
structure associated with the center-of-mass motion, but also their compos- 
ite nature in terms of the degrees of freedom of the constituent fermions. 

The key feature of ultracold Fermi atoms that has allowed the realization 
of the BCS-BEC crossover is the possibility of varying essentially at will the 
strength of the attractive interaction between fermions of different species 
[2j, attraction which results in the formation of Cooper pairs in a medium, 
on the one hand, and of composite bosons in vacuum, on the other hand, 
out of the two fermion species. [In the case of ultracold atoms, the spin of 
an electron is replaced by an analogous quantum number associated with 
the atomic hyperfine levels.] Owing to this unique possibility, ultracold 
atoms should be regarded as prototype systems, with respect to others in 
Nature for which this possibility is hindered. Specifically, in ultracold atomic 
gases the attractive interaction is varied through the use of the so-called 
Fano-Feshbach resonances, which are characterized by a resonant coupling 
between the scattering state of two atoms with near- zero energy and a bound 
state in a close channel [3]. Changing (through the variation of a static 
magnetic field B) the position of the bound state with respect to threshold 
in a suitable way, one can modify the value of the (fermionic) scattering 
length aF from negative values before the formation of the bound state in 
the two-body problem to positive values once the bound state is formed [I] . 
As an example, Figjl] shows the scattering length for the collision of two 6 Li 
atoms vs B. 

In this context, the dimensionless parameter (fepap) -1 acquires a spe- 
cial role for the corresponding many-body system at finite density. Here, 
the Fermi wave vector hp of the non-interacting system (which is defined 
as kp = yj2mEp both for a homogeneous and a trapped system - see be- 
low) is a measure of the (inverse of the) interparticle distance, m being 
the fermion mass and Ep the Fermi energy of non-interacting fermions (we 
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Figure 1: Scattering length for 6 Li atoms vs magnetic field. The inset 
amplifies the behavior of the narrow resonance. (Adapted from Fig. 2 of 
Ref.JS].) 

set % = 1 throughout). When the Fano-Feshbach resonance is sufficiently 
"broad" (like for 6 Li and 40 K atoms used in experiments thus far), in fact, 
the many-body fermion problem can be described in a simplified way by 
a single-channel Hamiltonian with an instantaneous short-range interaction 
[S]. The strength of this interaction, in turn, can be parametrized in terms of 
the above scattering length ap of the two-body (molecular) problem, which 
shares the same ultraviolet divergency associated with the short-range char- 
acter of the two-body potential O [7]. In this way one ends up with all 
physical quantities of interest for the many-body system depending on the 
interaction only through the parameter (kpap) . 

In terms of this parameter, one finds that for most physical quantities 
the crossover between the BCS and BEC regimes is exhausted, in practice, 
within a range ~ 1 about the unitary limit at {kpap) = where ap di- 
verges. Outside this limited range, the BCS and BEC regimes (whereby 
(A^of) -1 ~ —1 and 1 ~ (kpap)^ 1 , in the order) are characterized by 
the product kp\ap\ being quite smaller than unity (corresponding to a 
diluteness condition), so that theoretical approaches can in principle be 
controlled in terms of this small quantity in these two separate regimes. 
No such small parameter evidently exists, however, in the unitary regime 
about (kpap) -1 = 0, whose theoretical description consequently constitutes 
a formidable task. 
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It is then clear that theoretical treatments of the BCS-BEC crossover 
should provide as accurate as possible descriptions of the two regimes where 
the above diluteness condition applies, either in terms of the constituent 
fermions (BCS regime) or of the composite bosons (BEC regime). Specif- 
ically, this has to occur via a single fermionic theory that bridges across 
these two limiting representations, by recovering controlled approximations 
on both sides of the crossover and providing at the same time a continuous 
evolution between them, thereby spanning also the unitary regime where 
use of the theory could a priori not be justified. 

The prototype of this kind of approach is represented by the BCS the- 
ory itself at zero temperature. As remarked originally by Leggett (8] (see 
also Ref.[9j), the BCS wave function is quite more general than originally 
thought, in the sense that it contains as an appropriate limit the coherent 
state associated with a Bose-Einstein condensate of composite bosons made 
up of opposite-spin fermions. This limit is reached when the occupation 
numbers of all possible fermionic single-particle states are much less than 
unity, so that the Fermi surface is completely washed out. 

The argument can be made more quantitative by solving the coupled 
gap and density equations provided by the BCS theory [TOl [TJ] for a homo- 
geneous system at T = (for which an analytic solution exists in terms of 
elliptic integrals [12]), to obtain the gap (order) parameter Ao and fermionic 
chemical potential fiQ as functions of (A^cif) -1 as shown in Figj2j Note that 
the chemical potential crosses over as expected, from the value Ep of the 
Fermi energy of non-interacting fermions in the BCS limit, to (half the 
value of) the binding energy £o = {rna? F )~ l of the two-body (molecular) 
problem within the single-channel model in the BEC limit. In both limits, 
Aq/I/j-oI *C 1 albeit for different physical reasons. 
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Figure 2: (a) Gap parameter Ao and (b) chemical potential fj,o for a homo- 
geneous system at zero temperature vs the coupling parameter (/c^ap)" 1 , 
evaluated across the BCS-BEC crossover within mean field. 

The BCS theory is a mean-field approximation which relies on the Cooper 
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pairs being highly overlapping in real space [TO], so that their effects can be 
dealt with "on the average" . As such, it is expected to be a valid approxima- 
tion even at finite temperature whenever this condition is satisfied. It should 
accordingly apply to the BCS limit of the BCS-BEC crossover, but not to 
the unitary or BEC regimes where the typical length scale for correlation 
between two fermions with different spins becomes comparable with the in- 
terparticle spacing kp 1 . This is shown in Fig jj}] where the (zero-temperature) 
intra-pair coherence length £ pa ; r is plotted vs (kpap)^ 1 . 
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Figure 3: Pair coherence length £ pa i r (dashed line) and phase coherence 
(healing) lenght £ p hase (full line) vs (/c^of) -1 ; evaluated at zero tempera- 
ture as in Ref.[12] according to their definitions given in Refs.[T3] and |14j . 
respectively. 

It is then evident that, away from the BCS limit, inclusion of fluctuation 
corrections beyond mean field becomes essential to account for the relevant 
physical properties of the system. An equivalent way of stating the problem 
is that, away from the BCS limit, the intra- and inter-pair coherence lengths 
are expected to differ considerably from each other. This is also shown in 
Figj3] where the inter-pair coherence length £ p hase is reported for compar- 
ison. [In the BCS limit the two lenghts differ by an irrelevant numerical 
factor (£ pair ~ (3/\/2)£ P h ase ) owing to their independent definitions, so that 
the two curves in Figj3] are parallel to each other in this coupling regime.] 
In particular, in the BEC regime £ pa i r corresponds to the size of a compos- 
ite boson while £ pria se represents the healing length associated with spatial 
fluctuations of the center-of-mass wave function of composite bosons. 

It was indeed within this framework that the BCS-BEC crossover at- 
tracted attention also for high-temperature (cuprate) superconductors |T5l 
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[TB] . for which the product fcir£pair was estimated to be about 5-=-10 in con- 
trast with more conventional superconductors for which it is of the order 
10 3 -T- 10 5 . Several theoretical works were then put forward on the BCS-BEC 
crossover in this context [T71 [181 [131 032 GO], with the limitations, however, 
that the origin and characteristics of the attractive interaction at the basis 
of this crossover were not known for cuprate superconductors. These limita- 
tions have eventually been fully removed with the advent of ultracold Fermi 
atoms, to which we shall limit our considerations in the following. 

One related reason to invoke the inclusion of fluctuation corrections be- 
yond mean field stems from the values obtained within BCS theory for the 
critical temperature at which the order parameter vanishes. Only a numeri- 
cal solution of the coupled gap and density equations is amenable for generic 
values of (A^of) -1 , but analytic results can still be obtained in the BCS and 
BEC limits. One gets [6]: 

8 Ef e 1 

k B T c ~ ^- exp{vr/(2a F A; F )} (1) 

7r e 

in the BCS limit (where fee is Boltzmann constant and 7 Euler constant 
with e 7 /-7r ~ 0.567), and 

2 In (e /E F f 2 

in the BEC limit, respectively. While the result ([I]) corresponds to what is 
familiar from BCS theory for weak coupling [10] . the result ([2]) does not co- 
incide with what one would expect in the BEC limit, namely, the expression 

9 /q 

of the Bose-Einstein condensation temperature /cb^bec = 3.31 n B /(2m) 
where ub = n/2 is the density of composite bosons in terms of the density 
n of the constituent fermions. On the contrary, the expression ([2]) increases 
without bound when approaching the BEC limit for 1 << (kpap) ■ 

The points is that the critical temperature obtained from the solution 
of the mean-field equations corresponds to the process of pair formation 
and not of pair condensation. The two temperatures coincide only in the 
BCS (weak-coupling) limit |21j . where pairs form and condense at the same 
time. In the BEC (strong-coupling) limit, on the other hand, pairs form at 
a higher temperature than that at which they eventually condense owing to 
quantum effects. Accordingly, the expression ([2]) signals the phenomenon of 
pair dissociation, and as such it must be regarded as a "crossover" tempera- 
ture T* which does not correspond to a true phase transition. The complete 
plot of T* obtained by solving numerically the mean-field gap and density 
equations throughout the BCS-BEC crossover is shown for the homogeneous 
(h) and trapped (t) cases in FigjiJ where it corresponds to the upper dashed 
and full lines, respectively (the remaining two curves labeled by T c result 
instead beyond mean field and will be discussed in Section [2]) . 
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Figure 4: Temperature vs coupling diagram for the trapped (full lines) and 
homogeneous (dashed lines) system, where the critical temperature T c and 
pair-breaking temperature T* are shown. Each temperature is normalized 
to the respective Fermi temperature Tp. (Adapted from Fig.l of Ref.|22j.) 



It is thus evident from the above discussion that the main limitation of 
the mean-field description we have considered thus far is that it includes 
only the degrees of freedom internal to the pairs which are associated with 
pair-breaking, but omits completely the translational ones. The latter are 
responsible for the collective sound mode, which represents the main source 
of elementary excitations in the BEC regime [23J. To overcome this severe 
limitation for a sensible description of the BCS-BEC crossover in terms of 
a fermionic theory, it is then necessary to go beyond mean field and include 
pair-fluctuation effects as discussed in the next Section. 

2 Inclusion of pairing fluctuations 

A diagrammatic approach for fermionic pairing fluctuations was first consid- 
ered by Galitskii [24J for a dilute Fermi gas with strong short-range repulsion 
|25j . There it was shown that the relevant fermionic self-energy can be taken 
of the form depicted in FigjHJ where T° is the pair propagator describing the 
repeated scattering in the medium between two fermions of opposite spins. 

The short-range nature of the potential requires one to introduce at the 
outset a regularization procedure that eliminates the ultraviolet divergences. 
This is done by exploiting the two-fermion problem in vacuum, which shares 
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Figure 5: Single-particle fermionic self-energy in the normal phase (upper 
panel) expressed in terms of the pair (ladder) propagator T° between two 
fermions of opposite spins (lower panel). Full and dashed lines represent 
the fermionic propagator and interaction potential, respectively, while the 
labels k (A/) and q correspond to fermionic and bosonic four-vectors, in the 
order. 



the same sort of divergences and for which the (positive) strength vq of the 
repulsive interaction can be related to the ultraviolet cutoff k$ in wave- 
vector space through the following equation for the two-body t-matrix in 
the low-energy limit [26]: 

m 1 f ko dk m 
4^ = V Q + J (2tt)3 • (3) 

This relation defines the (fermionic) scattering length ap, which is positive 
in this case and remains smaller than the range %/{2ko) of the potential if 
its strength vq is kept finite. 

To the leading order in ap, this self-energy results in a repulsive "mean- 
field shift" (A-trap/m) n/2 of the chemical potential, where n is the total 
fermion density for both spin components. This is because, to the leading 
order in ap, T°(q) ~ (—Airap/m) and the loop in the upper panel of FigjS] 
gives half the fermionic density n (we consider throughout the case of equal 
populations of spin up and down fermions). Terms up to the second order 
in ap were also obtained by the Galitskii original approach [23]. 

The above choice of the self-energy emphasizes pair-fluctuation effects 
via the repeated scattering of two (opposite spin) fermions in the medium. 
As such, it has been considered physically relevant also to the case of an 
attractive short-range potential with a negative vq [TTJ [271 121] , for which 
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ap has the typical resonant behavior of Fig |T] associated with the BCS-BEC 
crossover. By this extrapolation, the formal structure of the Galitskii self- 
energy is carried over to the domain of strong coupling, and even further to 
the repulsive side of the resonance where a molecular state forms [29J. 

Let's write down explicitly the analytic expressions corresponding to 
the diagrams depicted in Fig(5j for the simplest case when all fermionic 
propagators appearing therein are "bare" ones |27| . One has: 

S(k,w n ) = -fc B T^ |^r (q,^)Go(q-k,^-L.; n ) (4) 

for the fermionic self-energy, and 
(—1) m f dk 



r°(q,n„) 4vra F J (2tt) 



k B T Y Go(k, w„)G (q - k, O 



111 

w 

(5) 

for the (inverse of the) pair propagator. Here, Go(k, uj n ) = [iu n — £(k)] 1 
is the bare fermion propagator (£(k) = k 2 /(2m) — fi being the free-particle 
dispersion measured with respect to the chemical potential fi), while oj n = 
Txk^T{2n + 1) (n integer) and Cl u = ^i^k^Tv [y integer) are fermionic and 
bosonic Matsubara frequencies at temperature T, in the order. Note how 
the strength vq of the attractive interparticle potential has been eliminated 
in the expression ^ in favor of the scattering length ap via the relation ([3]), 
which now admits also negative value for clf consistently with the behavior 
shown in Fig[T| 

With the self-energy Q one dresses the bare fermion propagator to ob- 
tain the full propagator 

G(k,w n ) = ^rj ^ — — r , (6) 

in terms of which the chemical potential can be eventually eliminated in 
favor of the density via the expression (rj = + ): 

r dk 

n = 2k B Tj2e Mnr >J p»G(k,, n ) . (7) 

On physical grounds, the relevance of the expressions Q and §5§ to the 
BCS-BEC crossover can be appreciated from the following considerations. 
While in the BCS weak-coupling limit (where ap < and (kpap) -1 <C 
— 1) the pair propagator maintains formally the same expression T°(q) ~ 
—Airap/m of the repulsive case, in the BEC strong-coupling limit (where 
< ap and 1 <C (kpap)' 1 ) it acquires the polar structure of a free-boson 
propagator [T51 [7] : 

r°(q,n v ) = — ( 8 ) 

m z ap i\l v — q z /(4m) + \ib 
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where the bosonic chemical potential [ib reduces to 2/x + Eq in this limit 
when the composite bosons have size ~ ap. In this limit, we may expand 
the fermionic propagator ^ to the lowest order in £ 

G(k,w n ) ~ G (k,o; n ) + G (k,o; n ) S(k,w n ) G (k,u; n ) , (9) 

and consistently approximate the self-energy Q in the form: 



E(k,w n ) ~ -G (-k,- Wn )A; B r^e in - |^r°(q,^) . (10) 
In this way, we obtain for the density ([7]): 



(2tt) ; 

dk 1 _ f dq , 



of f *1 _ ^--\ 

n 7 (2tt) 3 e «k)/(* B T) + 1 / (2ttF ^ 



(2vr) 3 e£W/(*BT) + i 7 (2vr) 3 " ^ ifi„ - q 2 /(4m) + /x B 

(11) 

where use has been made of the result 

/(^3^ T E G o(k,^) 2 G (-k,-o; n ) ~ (12) 

which is valid when [i ~ — £q/2 is the largest energy scale in the prob- 
lem. Under these circumstances, the first term on the right-hand side of 



Eq.(ll) is strongly suppressed by the smallness of the fugacity e^^ ksT \ 
while the second term therein represents the density ns of a non-interacting 
system of (composite) bosons with chemical potential /i B , yielding even- 
tually n ~ 2ub- With the inclusion of pairing fluctuations, the density 
equation ([7]) thus reproduces the standard result for the Bose-Einstein con- 
densation temperature ^bTbec = 3.31n^ 3 jms where tub = 2m is the mass 
of a composite boson. 



Note further that, if only the first term on the right-hand side of Eq.( 11 ) 
were retained, one would get for the chemical potential: 



n 

2 Im 



27T 

k B T) 



(13) 



which coincides with the classical expression at temperature T [25]. Setting 
in this expression fi ~ — £o/2 and T = T c , the value ^ for the critical 
temperature T c is readily recovered. 

Quite generally at any coupling across the BCS-BEC crossover, the crit- 
ical temperature is obtained from the normal phase by enforcing in Eq.([5]) 
the instability condition l/r°(q = 0, Q u = 0) = 0, in conjunction with the 
density equation ([7]). The resulting values for T c are plotted for the homo- 
geneous (h) and trapped (t) cases in FigJiJ where they correspond the lower 
dashed and full lines, respectively. In both cases, in the temperature window 
between T c and T* composite bosons are formed but not yet condensed. 
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In the context of the BCS-BEC crossover, pairing fluctuations in the 
normal phase were first considered by Nozieres and Schmitt-Rink (NSR) 
|30j . with the purpose of obtaining a sensible extrapolation of the critical 
temperature from the BCS to the BEC limit (in Ref. [SO] the density equation 
was obtained by an alternative procedure via the thermodynamic potential) . 
It was later remarked in Ref. [31 J that the NSR procedure corresponds to a 
t-matrix theory in which one keeps only the lowest-order terms of Eq. ^ for 
all couplings and not just in the BEC limit. In practice, differences between 
the numerical results, obtained alternatively by the NSR procedure or by the 
approach based on Eqs.Q-Q where the expansion ^ is avoided, remain 
sufficiently small even in the unitary region. 

The approach for the normal phase based on Eqs.Q-Q was consid- 
ered in Ref.[27J to study fermionic single-particle properties above T c in 
the homogeneous case, and later extended to consider the effects of a trap. 
Owing to the presence of two bare fermion propagator Go in the particle- 
particle bubble of Eq.([5]), this approach is sometimes referred to as the 
"Go-Go t-matrix". This is to distinguish it from alternative t-matrix ap- 
proaches, notably: (i) The "G-Go t-matrix" approach [32] where one bare 
Go and one self-consistent G enter the particle-particle bubble defining the 
pair propagator, while a bare Go is kept in the definition of the fermionic 
self-energy (cf. Eq.@); (ii) The "G-G t-matrix" approach [33j where all 
single-particle Green's functions are self- consistent ones. These alternative 
approaches were both utilized recently to study the fermionic single-particle 
spectral function in the normal phase [34|, 135] . [It should be mentioned in 
this context that a t-matrix approach formally similar to the Gq-Gq one was 
proposed in Ref. [36], where the bare value of the chemical potential for the 
non-interacting Fermi gas was inserted in the self-energy in the spirit of a 
1/N expansion.] 

While the G-Gq and G-G t-matrix approaches have been implemented 
according to their strict definitions, the pairing-fluctuation approach of Eqs.(Q 
([7]) can be allowed to retain the original flexibility of the diagrammatic 
fermionic structure which is "modular" in nature. In this sense, it can be 
progressively improved by including additional self-energy corrections which 
are regarded important, especially in the BCS and BEC regimes where the 
approximations can be controlled. This implies, in particular, that the pair 
propagator in the expression Q can be dressed via "bosonic" self-energy 
insertions, which lead, for instance, to the Gorkov and Melik-Barkudarov 
corrections [21] on the BCS side and to the Popov theory for composite 
bosons [37] on the BEC side. Consideration of the latter is expected to be 
especially important on physical grounds, since it effectively introduces a re- 
pulsive interaction among the composite bosons which ensures, in particular, 
the stability of the system under compression. 



10 




Figure 6: The isothermal compressibility dn/dfi (in units of n/Ep) vs T/Tp, 
as obtained at unitarity from the pairing-fluctuation approach based on 
Eqs.Q-Q, is shown to diverge at T c (full line). In contrast, the curve 
corresponding to the free Fermi gas goes smoothly through T c (dashed line). 



One major shortcoming of the pairing- fluctuation approach of Eqs.Q- 
is, in fact, that it leads to a diverging compressibility when the tempera- 
ture is lowered down to T c from the normal phase. This behavior is shown in 
FigJBJat unitarity, and can be ascribed to the fact that the pair propagator 
([5]) corresponds to non-interacting composite bosons. 

The price one has to pay, for setting up theoretical improvements over 
and above the pairing-fluctuation approach discussed in the present Section, 
is the unavoidable increase of their numerical complexity when calculating 
physical quantities. Some of these improvements will be discussed in the 
next Section. 

The above pairing-fluctuation approach can, in addition, be extended to 
the superfluid phase below T c , whereby the pair propagator acquires a matrix 
structure that maps onto the bosonic normal and anomalous propagators 
within the Bogoliubov theory |38} [39] . This extension (together with its 
Popov refinement [37]) will also be considered in the next Section. 
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3 Bogoliubov and Popov approaches, and the resid- 
ual boson-boson interaction 



A pairing-fluctuation approach was implemented on physical grounds be- 
low T c in Refs.[3Hl EH], by adopting a fermionic self-energy in the broken- 
symmetry phase that represents fermions coupled to superconducting fluc- 
tuations in weak coupling and to bosons described by the Bogoliubov theory 
in strong coupling. This approach has allowed for a systematic study of the 
BCS-BEC crossover in the temperature range < T < T c . 

A diagrammatic theory for the BCS-BEC crossover below T c was ac- 
tually first proposed by Haussmann [18 1, by extending the self-consistent 
t-matrix approximation to the broken-symmetry phase. While the ensuing 
coupled equations for the chemical potential and order parameter were ini- 
tially solved at T c only, an improved version of this self-consistent theory 
was recently implemented for the whole thermodynamics of the BCS-BEC 
crossover [33J. We postpone an explicit comparison with this alternative ap- 
proach to Section [4j where a selection of numerical results will be presented. 

By the approach of Refs.[38, 39], the pair propagator in the broken- 
symmetry phase has the following matrix structure: 




lute) r 12 ( g ) \ 

r 2 i(g) r 22 ( q ) J A(q)A(- q ) - B{ q f 



where 



At \ 171 f dk 



k B Tj2Sn(k + q)gn(-k) 



m 



(14) 



(15) 



iiraF J (27r) J 
r dk 

B(q) = Jj^k B Tj2Mk + q)G2l(-k). (16) 

This structure is represented diagrammatically in Figj7^i, where only com- 
binations with £l = H-'l and £r = £' R survive the regularization we have 
adopted for the potential [cf. Eq.([3])]. It represents an approximation to 
the Bethe-Salpeter equation for the fermionic two-particle Greens function 
in the particle-particle channel. In the above expressions, q = (q, Q u ) and 
k = (k, w n ) are four- vectors, and 



Gii{k,u) n ) = - — ^ = - g 2 2{-k, -U} n 



Qi 2 {k,u n ) = — 2 = ^2i(k,w n ) (17) 



are the BCS single-particle Green's functions in Nambu notation |llj . with 
E(k) = ^/£(k) 2 + A 2 for an isotropic (s-wave) order parameter A (which 
we take to be real without loss of generality) . 



12 



q-k q-k' q-k q-k' q-k q-k' 











r B "i > * > ; > 




T(q) 






+ 






in 







+ ... (a) 



k k' k k' k k' 




Figure 7: Single-particle fermionic self-energy for the broken-symmetry 
phase (panel (c)), expressed in terms of the pair propagator T with Nambu 
structure (panel (a)). The BCS contribution to the self-energy is shown in 
panel (d), and the corresponding self-energy for the normal-phase of Figj5] 
is also reported in panel (b) for comparison. (Reproduced from Fig.l of 
Ref.[39].) 



In analogy to what was done for obtaining the expression ([8]) in the 
strong-coupling limit, one can show that the pair propagator ( |14[ ) reduces 
in the same limit to the following expressions: 



rnGz) 



221, 



8vr n B + iVL v + q 2 /(4m) 



m 2 aF 



and 



where 



ri2(<?) = r 2 i(g) 



^s(q) 



8tt 



fJ-B 



m 2 a F E B (q) 2 - {i9. v f 



1 



+ Ms 



m| 



(18) 



(19) 



(20) 



V \2m B 

has the form of the Bogoliubov dispersion relation [25], Ms = A 2 /(4|//|) = 
2^u + £o being the corresponding value of the bosonic chemical potential. 
Apart from the overall factor —8ir /(to 2 of) (and a sign difference in the 
off-diagonal component [38]), the expressions (18) and (19) coincide, respec- 



tively, with the normal and anomalous non-condensate bosonic propagators 
within the Bogoliubov approximation 
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For any coupling, in Ref.[39] the corresponding fermionic self-energy was 
taken of the form: 

dq 



En (A;) = -E 22 (-fc) = -k B T^ j ^ Tii(q, fi„) gn(q - k, Q v - u) n ) 

Ei 2 (fc) = E 21 (fc) = -A (21) 

where En is shown diagrammatically in Figj7j3 (with £l = £' L = £r = £' R = 
1) and E i2 in FigjTfi. 

With this choice of the self-energy, the fermionic propagator is then 
obtained by solving Dyson's equation in matrix form: 



G u Hk) G^(k) 
G^(k) G^(k) 



G (k)- 1 




W E n (fc) Si 2 (fc) \ 
-Go(-k)- 1 J [ E 2 i(fc) S 22 (fc) J 

(22) 

Note that the BCS expressions (17) for the single-particle Greens functions 
result by neglecting in Eqs.(22) the diagonal self-energy terms associated 
with pairing fluctuations. 

By the approach of Ref . [39] . the "normal" propagator G\\ is inserted in 
the density equation 



f dk 

n = 2 k B T ^ e luJnV J —G 11 (k,co n ) 
which replaces Eq.Q below T c , while in the gap equation 



dk 

W) 



3 £l2(k,W n ) 



(23) 



(24) 



the BCS "anomalous" propagator (17) is maintained (albeit with modified 
numerical values of the chemical potential and order parameter that result 
from the simultaneous solution of Eqs.(23) and (|24[)). This ensures that the 
bosonic propagators (14) remain gapless. 

It is instructive to consider once more the BEC limit, whereby the diag- 
onal part of the self-energy acquires the following approximate form [39J: 

8vr 1 



En(k,w„) ~ — ^— - ^—jz-r n' B (T) . 

m z aF i0J n + ?(k) 



Here, 



n' B {T) 



ciq 



u|(q)6(^B(q))-u|(q)6(-£?B(q)) 



(25) 



(26) 



represents the bosonic noncondensate density, with the Bose distribution 
b(x) = {exp[x/(k B T)] — l} -1 and the standard bosonic factors of the Bo- 
goliubov transformation [25j: 

2 / \ 2 / \ 1 2m s 

Mq) = - 1 = — - 



+ fj, B ~ E B (q) 



2£ B (q) 



(27) 
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In this case, solution of the Dyson's equation (22) yields: 



Gn(k,u; n ) ~ (28) 



iu n - £(k) 



-«k) 



with the notation Ap = &im' B (T) / (m^ap). When inserted into the density 



equation (23) the above expression gives: 

n ~ ^ (A 2 + A 2 g ) = 2 (n (T) + n' B (T)) (29) 

where the condensate density uq(T) is identified via A 2 = 8imo(T)/(m 2 ap). 

It is relevant to comment at this point on the value of the scattering 
length as, which results in the BEC limit of the above approach from the 
residual interaction between composite bosons. This value is obtained, for 



instance, by manipulating the gap equation (24) in this limit, yielding: 
A 2 



2|/i|e - 2\fi\ \ ~ fi B • (30) 



With the relation between A 2 and uq utilized in Eq. ( 29 ) and the asymptotic 
result \fj,\ ~ (2ma F ) , the expression (30) can be cast in the form (j,b = 



47raBno/ m B that corresponds to the value of the Bogoliubov theory with 
as = lap. 

This result can also be interpreted diagrammatically as being associated 
with the lowest-order (Born approximation) value for the effective boson- 
boson interaction [7] . This is represented in Figj8^, and can be obtained from 
the following expression where all bosonic four- momenta qi (i = 1, •••,4) 
vanish: 

« 2 (0, 0,0,0) = k B Tj2 Jj^G (p) 2 Go(-p) 2 



/ 9 \ 2 

dp 1 I map\ ( \'Kap\ 



(2vr) 3 4£(p) 3 I 8vr I \ m 



) (31) 



with the last line holding in the BEC limit. Apart from the overall factor 
{m 2 aF /(87r)) 2 (that compensates for the presence of the factor — 8ir/ (m 2 ap) 



in the expression (|8j) of the free-boson propagator), the result (31 ) is indeed 
consistent with a residual bosonic interaction corresponding to as = lap. 

The correct value for as (= 0.6of), which includes all possible scat- 
tering processes between two composite bosons in isolation, was originally 
determined in Ref.[40j from the exact solution of the Schrodinger equation 
for dimer-dimer elastic scattering, and later confirmed in Ref.jH] through 
a completely diagrammatic treatment at zero density. In this context, even 
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p+q 4 




(a) 



(b) 




+ - 



Figure 8: (a) Effective boson-boson interaction u-i. (b) Additional terms 
associated with the t-matrix i B for composite bosons. Light lines stand for 
free-fermion propagator and broken lines for fermionic interaction potential. 
Spin labels are not shown explicitly. (Reproduced from Fig. 2 of Ref. |37j ■) 



before this exact result was available, it was shown in Ref. [7] that the scatter- 
ing processes corresponding to the t-matrix diagrams for composite bosons 
(the lowest ones of which are depicted in Figj8]3) lead by themselves to 
a considerable reduction of the value of a B (~ 0.75ap) starting from the 
value as = 2clf of the Born approximation. As a matter of fact, the 
complete diagrammatic treatment of Ref. [41 J (that yields the exact value 
as = O.Gap) adds to the diagram of Fig|8^ all other additional (zero den- 
sity) processes which are irreducible with respect to the propagation of two 
composite bosons, and then uses the result in the place of the diagram of 
Fig|8^ as the new kernel of the integral equation depicted in FigJHJa. 

The above considerations suggest us a way to improve on the Bogoliubov 
approximation for composite bosons, in order to include the diagrammatic 
contributions leading to a refined value of a B with respect to the Born 
approximation (in the following, we shall limit ourselves to recovering the 
value a B — 0.75clf in the BEC limit). To this end, we first approximately 
obtain the pair propagators Tb for any value of the fermionic coupling, by 
adopting the following Dyson's type equation in matrix form |37j in the 



place of the expressions (14): 



T B (q) = T° B (q) + T° B (q)^ B (q)T B (q) . (32) 
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Here, Tg(q) is the free-boson propagator with inverse 



where r°(g) _1 is given by Eq.^, and 

v _ 42 / -2u 2 (0,g,0,g) iZ 2 (0,0,-g,g) \ , . 

Sb(9) " A I, U2(0,0,- g , ff ) -2u 2 (0,g,0,g) ) (34) 

is the bosonic self-energy within the Bogoliubov approximation, which con- 
tains two degenerate forms of the effective boson-boson interaction [cf. Figj8^]: 

«2(gi,g2,<?3,94) = k B TY^ Jj^^Go(-p)Go(p+q2)G (-p+qi-q4)G (p+q i ). 

(35) 

To guarantee the ladder propagators Ts(q) of Eq.(32) to be gapless when 
q = for any value of the fermionic coupling, we impose the condition: 

T ^ = 0)- 1 - = 0) - T}i{q = 0) = (36) 

which plays the role in the present context of the Hugenholtz-Pines theorem 
for point-like bosons |25j . The Bogoliubov approximation for the composite 
bosons in the BEC limit with an improved value of as then results [7], 
by replacing in Eq.(34) the boson-boson interaction u 2 with the following 
expression of the t-matrix for composite bosons (cf. Fig{8| 



(27^ (37) 

xu 2 {qi,q 2 ,q5,qi + <?2 - g5)F (<7 5 )r (gi + q 2 - q->)i B {qi + <?2 - <?5, 95, <?3, 94)- 

Further improvements can be implemented by replacing the pair propa- 
gators (14) with more refined descriptions of composite bosons in the BEC 
limit, and then using these improved descriptions throughout the BCS-BEC 
crossover to modify the fermionic single-particle self-energy accordingly. An 
example is the so-called Popov approximation for composite bosons |37| . 
whereby the bosonic self-energy represented diagrammatically in Figj9] is 
employed to modify the original Bogoliubov propagators ( 14 ). In the broken- 
symmetry phase, this bosonic self-energy has the form: 

x G U (P + q? Gn(-P) G 11 (q' -q-p) r n (q') . (38) 
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Figure 9: Graphical representation of the Popov self-energy for composite 
bosons, that results upon dressing the upper fermionic line in the particle- 
particle channel. An analogous dressing done for the lower fermionic line ac- 
counts for the factor of two in Eq,(p8|). (Reproduced from Fig. 4 of Ref. [37J.) 



The Popov propagators for composite bosons are then obtained as fol- 
lows, in terms of the corresponding Bogoliubov propagators (|14|): 



rfr p (g) rf 2 op ( q ) 



A(-q) 



- ^ op ( 

B(q) 



111 



B(q) 

v Pop 



n 



r% p (q) r 2 T(g) 



Pop/ 



[A(q) 



^Pop 



(q)n][A(-q) 



-q) 



ii 



B(q? 
(39) 

where A(q) and B(q) are given by Eqs.(15) and (16), in the order. The 
propagators (39) are gapless provided 



A(q = 0) - Y, P B op (q = 0)n - B(q = 0) = . 



(40) 



This generalizes to the present context the condition A(q = 0)—B(q = 0) = 
for gapless Bogoliubov propagators, and effectively replaces the gap equation 
(24) for all practical purposes. 

In addition, the same treatment that was made above to improve on the 
relation as = Zap in the BEC limit can be applied here, by first rewriting 
( 38 ) in terms of the bare boson-boson interaction ( 35 ) 

(g)n ~ - 2 k B T V f ^ u 2 (q\ q, q', q) T n (q') , (41) 



^Pop 



and then replacing u 2 by the t-matrix ts for composite bosons of Eq.(37). 
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The final form of the fermionic self-energy is eventually obtained by 
reconsidering the expressions (21), where now the Popov propagator r-Q^ 5 
of Eq.(39) takes the place of Tu while A satisfies the condition (40) in the 
place of the original gap equation. 

From a physical point of view, the relevance of the Popov approximation 
results because it introduces an effective repulsion among the composite 
bosons through the presence of their noncondensate density. The importance 
of this repulsion should be especially evident in the normal phase, when 
the Bogoliubov propagators (14) reduce to free-boson propagator ^ and 
miss accordingly this residual bosonic interaction. While commenting on 
FigJB] we have already pointed out that this is the reason for a diverging 
compressibility at T c when only "bare" pairing fluctuations are considered. 

In the next Section we shall discuss a number of thermodynamic as well 
as dynamical results obtained by implementing the Popov approximation 
in the normal phase throughout the BCS-BEC crossover (the unitary limit 
will specifically be considered). In this case, the (inverse of the) Popov 
propagator for composite bosons is obtained from the relation r Pop (g) _1 = 
r°(g) _1 — Y; B ° p (q), where S^ op is given by Eq.(41) with T replacing Tu. 
In addition, to improve on the description of the boson-boson scattering, 
we shall replace the bare v,2 in Eq.(41) by the t-matrix is for composite 
bosons given by Eq.(37). In this way, the 2-boson scattering will be dealt 
with beyond the Born approximation. 

In this context, it will be relevant to compare the results obtained by 
the above approach in the normal phase for thermodynamic and dynamical 
quantities |42j . with those obtained by an alternative approach based on a 
self-consistent t-matrix approximation, as described in Ref.[33j for the ther- 
modynamics and in Ref. [35] for the dynamics of the BCS-BEC crossover, 
respectively. [Results for a homogeneous system will only be presented.] In- 
terest in this comparison is also justified on physical grounds, by considering 
the different treatments of the effective boson-boson interaction which result 
from the two approaches. As remarked already, the Popov approach with 
£b replacing U2 concentrates on 2-boson scattering beyond the lowest order 
(Born) approximation, while the self-consistent t-matrix approach includes 
a sequence of 3, 4, • • •, n-boson scattering processes, where each process is 
dealt with at the lowest order. Although these alternative sets of processes 
(namely, improved 2-boson vs n-boson scattering) can be clearly identified 
by a diagrammatic analysis in the (BEC) strong-coupling limit [7], the ques- 
tion of how the relevance of these processes extends to the unitarity limit 
remains open and can be addressed only via numerical calculations. This 
question will be partially addressed in the next Section. 
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4 Results for thermodynamic and dynamical quan- 
tities 



Physical quantities that can be considered for a quantum many-body sys- 
tem are conveniently organized as single- and two-particle properties, and 
are correspondingly obtained in terms of single- and two-particles Green's 
functions. In addition, these properties may refer to the equilibrium state of 
the system or to excitations over and above this state. In the first case they 
can be conveniently obtained within the Matsubara formalism with discrete 
imaginary frequencies, while in the second case a (sometimes nontrivial) an- 
alytic continuation to the real frequency axis is required [25] . In the present 
context of a pairing-fluctuation diagrammatic approach to the BCS-BEC 
crossover, we shall limit ourselves to considering the chemical potential and 
the total energy per particle as examples of thermodynamic properties, and 
the single-particle spectral function as an example of dynamical properties, 
for which consideration of pairing fluctuations appears especially relevant. 

This relevance is most evident in the normal phase, because the occur- 
rence of pairing fluctuations acts to extend above T c characteristic effects 
of pairing (notably, what is referred to as the "pseudogap physics" associ- 
ated with the noncondensate density like in Eq.(|28|)), effects which would 
otherwise be peculiar of the broken-symmetry phase below T c only. 



4.1 Thermodynamic properties 

For a homogeneous Fermi gas in the normal phase, the fermionic chemical 
potential \i can be obtained from the density equation ([7]) and the total 
energy per particle from the following expression |25j : 

- = -k B T £ e-'/^ ^ + , + «J G(k,„„) (42) 



where N is the total particle number. In Eqs.Q and (42), different ap- 
proximations are embodied in different forms of the fermionic single-particle 
Green's function G. 

In particular, we shall consider approximate forms of G obtained within: 
(i) The t-matrix approach given by Eqs.Q-Q; (ii) Its further simplification 
(sometimes referred to as the Nozieres-Schmitt-Rink (NSR) approximation) 
whereby the fermionic propagator is expanded like in Eq. ([9]) for any coupling 
(and not just in the BEC limit); (hi) The Popov approach with an improved 
description of the boson-boson scattering as discussed in Section |3j (iv) The 
fully self-consistent (sc) t-matrix approach of Ref.[33] that was mentioned 
in Section [2] 

At unitarity, the results of these diagrammatic approaches for \x and 
E/N can also be compared with Quantum Monte Carlo (QMC) calcula- 
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Figure 10: (a) Chemical potential and (b) energy per particle (in units of 
Ep) vs the temperature (in units of Tp), as obtained: In Ref. [42] by the 
t-matrix approach (dashed lines), the NSR approximation (dashed-dotted 
lines) , and the Popov approach (thick full lines) ; In Ref . [33] by the fully self- 
consistent t-matrix approach (dotted lines); In Ref. |45j (full squares) and in 
Ref. [46] (full circles) by QMC calculations. Results obtained by the modified 
virial expansion of Ref. [47] (thin full lines) are also shown for comparison. 



tions, which are available over a wide temperature range. This comparison 



is shown in Fig. 10 Several features are here apparent. At high enough 
temperatures, all data progressively merge to the t-matrix approach, which 
is known to become exact in this limit where it reduces to the virial ex- 
pansion of Beth and Uhlenbeck [43]. While only minor differences appear 
between the t-matrix and the NSR approaches (with independent NRS cal- 
culations yielding comparable results [44J), the Popov approach is seen to 
a add positive contribution both to /i and E/N. This is in line with the 
expectation that the Popov approach takes into account the residual (re- 
pulsive) interaction among composite bosons |37j . which is missed by the 
t-matrix approach. In addition, the fully self-consistent t-matrix approach, 
which at high enough temperature should also asymptotically reduce to the 
t-matrix approach without self-consistency, shows deviations from the t- 
matrix and Popov approaches that are more marked in \i than in E/N . No 
compelling conclusions can, therefore, apparently be drawn by comparing 
self-consistent vs non-self-consistent pairing-fluctuation approaches as far as 
the thermodynamic quantities are concerned. Quite generally, it can be 
stated that good overall agreement results by comparing QMC calculations 
with diagrammatic pairing-fluctuation approaches, signifying that the latter 



are able to capture the relevant physical processes. Note finally from Fig. 10 
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t-matrix 


NSR 


Popov 


sc t-matrix 


QMC I 


QMC II 


p 


-0.57 g2] 


-0.599 [H] 




-0.64 [33] 


-0.63(5) gH] 




M (T = 0) 


0.455 !49 


0.401 [44] 




0.36 [33] 


~ 0.42 [45] 




T c 


0.243 |42| 


0.225 [42] 


0.199 [42] 


0.160 [33] 


< 0.15 [48] 


0.152(7) [46] 




0.365 [42 


0.336 [H] 


0.464 42 


0.394 [33J 


0.43(1) 48 


0.493(14) 46 


E(T C )/N 


0.296 [42] 


0.278 [42] 


0.337 [42] 


0.304 [33] 


0.270(6) 48 


0.31(1) [46] 



Table 1: Thermodynamic quantities obtained theoretically for a homo- 
geneous Fermi system at unitarity by alternative pairing-fluctuation ap- 
proaches and QMC calculations (relevant references are specified). Energies 
are in units of E F and temperatures of Tp. Here, the dimensionless param- 
eter P results from the relation 1 + /3 = E(T = 0)/E n [(T = 0) that holds 
at unitarity, where E n i is the energy of the corresponding non-interacting 
system. 



that the progressively increasing differences between the t-matrix and Popov 
approaches when lowering the temperature reflects the fact that the diver- 
gence of the compressibility resulting from the t-matrix approach (shown in 
Figj6]) is suitably cut off by the Popov approach, which yields a finite value 
for this quantity at T c [32]. 

Alternative theoretical approaches yield different values of the critical 
temperature T c , as shown in Table [T] These values can be compared with 
the corresponding ones that are extracted from experiments, as reported in 
Table [2] Although for this quantity the self-consistent t-matrix approach 
seems to perform better than the non-self-consistent one(s), one should be 
aware of the fact that additional corrections to the pair propagator r°, like 
those introduced in Ref . [21 | in the weak-coupling (BCS) limit to represent 
the medium polarization and shown to have a sizable effect on T c in that 
limit, might still act to reduce somewhat further the value of T c even at 
unitarity. Definite comparison with the experimental values of T c reported 
in Table [2] should then await for a proper inclusion of these additional cor- 
rections. To elicit a more quantitative comparison among theoretical and 
experimental thermodynamic quantities, Tables [T] and [2] list, in addition, the 
values of /i and E/N that are available both at zero temperature and T c . 

The sizable effects that pairing fluctuations have on the thermodynamic 
quantities of Fig 10 over and above the free- Fermi gas behavior can be ap- 
preciated by sketching therein the plots of [i and E/N for the non-interacting 
Fermi gas (recall, in particular, that fi n i(T = Q)/Ep = 1, E n \(T = Q)/(NE F ) - 
0.6, /i ni (T = 0.6T F )/E F = 0.625, and E ni (T = 0.6T F )/(NE F ) = 1.15). The 
effects of pairing fluctuations in these quantities are thus seen to extend 
over a wide temperature range up to several times T F , being related to the 
high-energy scale introduced in Ref. [55] in terms of the trace of the pair 
propagator T . 
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Exp. [50] 


Exp. [H] 


Exp. [52] 


Exp. 07] 


Exp. [53] 


Exp. [54] 


p 


-0.62(2) 


-o.68±g:iS 


-0.54™ 


-0.58(1) 


-0.54(5) 




M (T = 0) 


0.38(2) 


o.32±g:ig 


0.46™ 


0.42(1) 


0.46(5) 




T c 








0.157(15) 




- 0.15 










0.49(2) 







Table 2: Thermodynamic quantities obtained experimentally for a homoge- 
neous Fermi system at unitarity (references are specified). Energies are in 
units of Ep and temperatures of Tp. The chemical potential at T = is 
here obtained via the relation fi(T = 0) /Ep = 1 + j3 that holds at unitarity, 
the parameter (3 being directly measured. [Experimental data for E(T C )/N 
are not available for comparison with the theoretical values of Table [T]] 



There exists, however, an additional energy scale (usually referred to 
as the pseudogap) which is also related to pairing fluctuations but is in- 
stead characteristic of the low-energy physics about T c . This energy scale is 
most evident when looking at the properties of the single-particle spectral 
function, to be considered next. 

4.2 Dynamical properties 

The spectral function A(k, to) for single-particle fermionic excitations re- 
sults after analytic continuation of the fermion propagator G(k,u n ) from 
the Matsubara (io n ) to the real (to) frequency axis, via the relation A(k, to) = 
-Im G R (k, to) /vr where G R (k, to) is the retarded fermion propagator. Through 
a related analytic continuation of the fermionic self-energy S, A(k, to) can 
be eventually cast in the form: 

A k, W = — — \ ' > 43 

^ — ?k — Rc2 J (k, oo)\ z + [ImL(k, ui)\ z 

where again ^(k) = k 2 / (2m) — /x. For any given wave vector k, the frequency 
structure of the real and imaginary parts of S(k, to) thus determines the 
positions and widths of the peaks in A(\z,co). 

The archetype of a pairing-gap behavior for A(k, to) is embodied in the 



two-peak structure of the following expression (cf. Eqs.(17)): 



A(k, to) = u(k) 2 5(to - E(k)) + v(k) 2 6(to + E(k)) (44) 

where E(k) = ^C(k) 2 + A 2 and v(k) 2 = 1 - n(k) 2 = (1 - f(k)/£(k))/2, 
which holds at the mean-field level in the broken-symmetry phase. When 
pairing fluctuations beyond mean field are included [27], a two-peak struc- 
ture still persists in the normal phase above T c , although with broad and 



asymmetric peaks replacing the delta spikes of Eq. ( 44 ) while the total area 
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Figure 11: Peak positions of the spectral function at T c for negative (lower 
branch) and positive (upper branch) energies versus the wave vector when 
(fcFair) -1 = —0.72. The spectral function is here obtained within the t- 
matrix approach of Ref. |27j. Full lines represent a BCS-like fit. (Adapted 
from Fig. 15 of Ref. [27].) 



remains unity. Even in this case, the positions of the two peaks in the spec- 



tral function follow quite closely the BCS-like dispersions ±y£(k) 2 + A^ g , 
provided the value A pg of the pseudogap is inserted in the place of the BCS 
gap A of Eq.(44). An example of this behavior is shown in Fig 11 for weak 
coupling. 

A systematic study of the single-particle spectral function in the normal 
phase across the BCS-BEC crossover was originally performed in Ref. |27j 
within the t-matrix approach given by Eqs.Q-Q. Interest in this study was 
recently revived by the advent of a novel experimental technique for ultracold 
Fermi gases |56] , whereby the wave vector of photo-excited atoms is resolved 
in radio-frequency spectra taken at different couplings and temperatures. 
One should mention in this context the comparison made in Ref. |34J between 
theoretical results obtained by the t-matrix and G-Gq approaches, as well as 
the calculation performed in Ref. [35] within the fully self-consistent t-matrix 
approach. 

Similarly to what was done in subsection |4.1| for thermodynamic prop- 
erties, here we compare the results for A(k, oo) obtained alternatively by the 
t-matrix, the Popov, and the fully self-consistent (sc) t-matrix approaches 
(while referring to Ref. |42j for a more complete analysis of this compari- 
son). We then show in Fig. 12 the results obtained for A(k,co) at unitarity 
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Figure 12: Single-particle spectral function at unitarity when T = T c and 
k = kp, as obtained by: (a) The t-matrix approach of Ref.[27J; (b) The 
Popov approach of Ref.[32]; (c) The sc t-matrix approach of Ref.[35]. (The 
plot of panel (c) has been extracted from Fig.4 of Ref. [35J.) 



25 



and k = kp by the three approaches, at the respective values of the critical 
temperature. [Analytic continuation from Matsubara to real frequencies has 
been performed in panel (a) by the direct substitution iuj n —> oj+irj, in panel 
(b) by the Pade approximants, and in panel (c) by the maximum-entropy 
method.] Note how the two-peak structure that is evident in panel (a) re- 
mains noticeable in panel (b), but has essentially disappeared in panel (c). 
This is consistent with a general understanding [57] that non-self-consistent 
calculations favor pseudogap behavior while self-consistent calculations tend 
to suppress it. As the other side of the medal, one would tend to attribute 
|35j to self-consistent calculations a more precise description of thermody- 
namic properties with respect to non-self-consistent approaches. 
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Figure 13: Dispersion of the two peaks of the spectral function at unitarity 
and T = T c obtained by the t-matrix approach (empty squares), the Popov 
approach (full dots), and the sc t-matrix approach (stars). (Stars have been 
extracted from the curves in Fig. 4 of Ref.[35].) 



In this respect, it should be noticed that the Popov appproach, while 
improving considerably on the thermodynamic description with respect to 
the t-matrix approach as discussed in subsection |4,1[ preserves also an ev- 
ident pseudogap behavior in the single-particle spectral function. This is 
made evident in Fig, 13, where the dispersion of the two peaks of A(k.,u) is 
reported at unitarity and T = T c for the t-matrix, the Popov, and the sc 
t-matrix approaches. (Close to &f, where the two peaks in A(k, u) are broad 
and overlap each other, the position of the less intense peak was determined 
by subtracting from ^4(k, to) the profile of the most intense peak.) The value 
of the pseudogap, identified by (half) the minimum energy separation be- 
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tween the upper and lower branches, remains essentially unmodified when 
adding the Popov on top of t-matrix fluctuations, but it closes up when full 
self-consistency is included. 

Stringent comparison with both experimental data and QMC calcula- 
tions will eventually decide what version of pairing-fluctuation theories is 
able to provide the closest agreement for thermodynamic as well as for dy- 
namical quantities, in systems like ultracold Fermi gases where only the 
mutual attractive interaction can be responsible for their physical behavior. 

5 Concluding remarks 

A gas of ultracold Fermi atoms, whose mutual interaction is governed by a 
(broad) Fano-Feshbach resonance, represents a physical system of fermions 
containing only pairing degrees of freedom. This feature naturally conveys 
their theoretical description in terms of "pairing fluctuations" of several 
kinds in the particle-particle channel, which extend characteristic two-body 
processes to a finite-density situation. The difficulty here is that, at fi- 
nite density, the relevant processes of the pairing type can be unambigu- 
ously identified only in the weak-coupling (BCS) regime where a fermion 
description is appropriate and in the strong-coupling (BEC) regime where 
a description in terms of composite bosons holds, because in both regimes 
the presence of the small parameter fcp|ai?| guides the selection of the dia- 
grammatic contributions for dilute systems. In additions, in these regimes 
useful analytic approximations can be quite generally derived from these di- 
agrammatic contributions, which help considerably one's physical intuition 
in picturing the involved processes. This kind of physical intuition is hard to 
emerge from more numerically oriented approaches (like QMC calculations) 
or more abstract approaches (like the renormalization group methods [58j). 
thus making diagrammatic approaches to the BCS-BEC crossover more ap- 
pealing in this respect. In principle, diagrammatic approaches may not be 
controlled in the unitary region, which is intermediate between the BCS and 
BEC regimes and where the diluteness condition does not apply owing to 
the divergence of |of|. Nevertheless, the good control of the approximations 
which can be achieved separately in the BCS and BEC regimes and the lim- 
ited extension of the unitary region (—1 ~ (A^of) -1 ~ +1) enable one to 
formulate a single fermionic theory that bridges across the BCS and BEC 
regimes and is able to furnish a good description of the unitary region for 
most practical purposes. This is the spirit with which diagrammatic pairing- 
fluctuation approaches to the BCS-BEC crossover have been formulated and 
applied to a variety of problems, both in the normal and superfluid phases. 

It is finally relevant to mention that the interest in the physics brought 
about by consideration of pairing fluctuations is not limited to a system of 
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ultracold Fermi atoms. In particular, the issue of the possible occurrence 
of a pseudogap in single-particle excitations is of considerable interest both 
in condensed matter [59J and nuclear physics |60j . This gives to ultracold 
Fermi gases the role of prototype systems, in which issues of general interest 
can be conveniently addressed by exploiting the unprecedented flexibility 
that they provide in the control of their physical parameters. 
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